Cantori and dynamical localization in the Bunimovich Stadium 
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Classical and quantum properties of the Bunimovich stadium in the diffusive regime are reviewed. 
In particular, the quantum properties are directly investigated using an approximate quantum map. 
Different localized regimes are found, namely, perturbative, quasi-integrable (due to classical Can- 
tori), dynamical and ergodic. 



PACS numbers: 05.45+b 

I. INTRODUCTION 

Long time ago, some physicists in the Siberian winter 
were playing dice with chaos and quantum mechanics. 
This game resulted in & paper which has been quoted 
N+l times in literature^ and which is known as a mile- 
stone in Quantum Chaos. 

Far from being an isolated branch of physics, Quan- 
tum Chaos has revealed its importance in the last twenty 
years in many important physical applications of differ- 
ent fields: Solid State, Nuclear, Atomic and Mesoscopic 
physics, just to give few examples. 

The common paradigma in Quantum Chaos is the so- 
called Kicked Rotator Model (KRM) whose dynamics is 
described by the Chirikov standard map (CSM). It rep- 
resents the safe retreat of many researchers working in 
Quantum Chaos. Indeed, even if still now papers dealing 
with some hidden property of the Chirikov Standard Map 
appearcl, (after the monumental work of Boris Chirikov 
published in 19790), or regarding new mathematical ad- 
vances, in the knowledge of the Floquet spectrum of the 
KRMcl, we may say that its general behavior is quite well 
understood, at least from the physical point of view. 

In this paper we show how a different physical problem 
can be explained using old results borrowed from KRM. 

The model under current investigation is the Buni- 
movich StadiumO. The properties of this model are well 
known in literature, both from physical and mathemati- 
cal point of view. Here we are interested in this particular 
billiard, characterized by a straight line 2a much smaller 
than the semicircle radius R : e = a/R -C 1. Preliminary 
studiesO have shown that the classical motion in the an- 
gular momentum is diffusive and it can be conveniently 
described by a 2-d area preserving map. In the same 
papertl the study of the nearest neighbor level spacing 
distribution (NNLSD) exhibits a different behavior de- 
pending on the energy range where the statistics is taken. 
In particular, at fixed eC 1 the NNLSD shows a smooth 
crossover from a Poisson to a Wigner-Dyson distribution 
as the energy is increased. The borders below and above 
which one can expect a particular distribution were thp 
oretically predicted and confirmed by numerical data 



The region characterized by intermediate statistics^ 
be associated, on the basis of a well defined pictureE 
with the presence of dynamical localization. In a sense, 
the qualified adjective "dynamical" was, at that stage 
quite inappropriate. Indeed only and indirect proof of 
dynamical localization were given, based on level statis- 
tics. The first example of localized eisenstates, in the 
angular momentum basis, was given inu for a rough de- 
formation of a circular billiard. In this case, due to the 
finite number of harmonics describing the smooth mod- 
ification of the boundary, direct exponential localization 
was found and the equality between quantum localization 
length and classical diffusion rate established. 

In this paper we enforce this viewpoint by studying 
directly the quantum dynamics instead of eigenfunctions 
and eigenvalues. This can be efficiently realized, from 
the numerical point of view, only quantizing the classi- 
cal map. The obtained numerical datal!3 indicate that 
quantum equilibrium distributions, differently from the 
classical ones, are algebraically localized in the angular 
momentum space. Further numerical investigational^ of 
the stadium eigenfunctions confirmed their algebraic lo- 
calization in the in angular momentum space. Neverthe- 
less it is possible to distinguish among different quantum 
regimes. 

Besides the perturbative regime, characterized by triv- 
ial periodic dynamical behavior, other different kind of 
localization can be identified. The first one, called quasi- 
integrable, has been associated to classical Cantori. An- 
other one is marked by dynamical localization, despite 
the algebraic localized distribution (with this word we 
mean a situation in which classical diffusion rate and 
quantum localization length hav|e-,the same numerical 
value). Recent analytical studiestU support this analy- 
sis. 

The paper is organized as the following: in section II 
we consider the classical dynamics inside the Bunimovich 
stadium as given by a suitable perturbed twist map. 
Classical properties of the discontinuous map are then 
investigated in section III, while section IV is devoted 
to its quantum properties. Finally in section V the es- 
timates about different borders in the energy-parameter 
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plane are summarized. 



with the Chirikov Standard Map, for which J — J ~ 
K sin 9 is a continuous function in the interval [0, 2it). 



II. MAPPING THE BUNIMOVICH STADIUM 



The ensemble dynamics of classical point particles hav- 
ing energy E, unit mass and initial angular momentum 
Iq, colliding elastically inside the Bunimovich Stadium, is 
described, when e <C 1, by the following map (R — Ip 



I 



Z - 2e sin 6> sgn(cos 9)^2E~l% 

9 + n - 2 asm(l/V2E) mod-27r 



(1) 



In (0) 9 is the angle measured from the center of the 
Stadium, I the angular momentum measured from the 
same center and the overlined variables (#, I) indicate the 
values taken after the collision with border. This map has 
been obtained by neglecting collisions with the straight 
lines, terms 0(e 2 ) and in the local approximation (small 
variations in the angular momentum). Indeed, it is easy 
to understand that while in the first equation of ([!]) I can 
grow infinitely, the second equation loses its validity when 
\l\ approaches its maximum value \/2E. The existence 
of such a bound is a trivial consequence of the energy 
conservation. 

If we now put Iq = and approximate asin(ir) with its 
argument (since we must exclude the values \x\ ~ 1), we 
get: 



I — I + k sin 9 sgn(cos 9) 
9 = 9 + Tl mod-27r 



(2) 



where we put k = 2e\/2E, T = 2/V2E, I -> -I and 
7r has been neglected since-,sin0 sgn(cos#) is 7r-periodic. 
Map (]]]) has been obtained^ for e«l, therefore (||) holds 
when kT = 4e < 1. The case kT > 1, which represents a 
possible regime for (||) has no physical meaning here and 
it will not be taken into account. 

Map (||) on the cylinder [0, 2tt) x (— oo, oo) has been 
recently investigated!!!] for kT < 1. Results can be then 
extrapolated to our case assuming \l\ <C \J2E. In the 
next sections we present a detailed study of classical and 
quantum properties of map (0). 



III. DISCONTINUOUS TWIST MAPS 

Let us write (^) by introducing the variables J = IT, 
K = kT in the following way : 

J = J + K sin 9 sgn(cos 9) 

9 = 9 + J mod-27r (3) 

such that we have single parameter K. Map (||) be- 
longs to a particular class of discontinuous twist maps. 
We use the word "discontinuous" to mark the difference 
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FIG. 1. Poincare surface of section for map (3) and 
K = 0.001. One particle has been iterated 2 • 10 4 times. 

The most studied case in the set of discontinuous func- 
tions is the saw-tooth map (STM) where the change in 
the angular momentum is given by J — J = K (9 — ir)/ir. 
For this map the classical transport properties have been 
studied, quite long agoO. The relevant difference of the 
STM, if compared with the CSM, is characterized by the 
absence of a KAM structure. Indeed, since the hypothe- 
sis of KAM theorem are not satisfied, we do not expect 
KAM tori exist for any K value, independently from its 
smallness. Namely, differently from the CSM, where for 
K < 1 the motion is typically regular on invariant tori, 
here one finds absence of KAM tori for any K; neverthe- 
less, invariant structures still exist. Indeed Cantori can 
be defined in the same way as in the STMLL3. An example 
of the motion in the neighborhood of Cantori is shown 
in Fig.^. As one can see the motion is far from being 
chaotic even if, upon the increasing of time, a single or- 
bit can explore, in a dense way, the whole phase space. 
Moreover, the resulting motion is conveniently described 
by a diffusive equation for the distribution function as 
can be inferred by looking at the behavior in time of the 
average squared momentum (see Fig.||a). 

The linear growth in time, after a transient time (see 
Fig.^) and the corresponding Gaussian distribution in 
angular momentum, at a given time, (see for instanceO) 
are usually taken as a. common reference for the existence 
of a diffusive motiond. 
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FIG. 2. Growth of the average momentum spreading in 
time for the discontinuous map (3) at K = 0.001. b) is a 
magnification of the small area indicated among dashed lines 
in the left corner of a). The initial ensemble consists of 2000 
particles having the same momentum Jq = 0.1 and random 
phases 9. 

For such discontinuous maps, the resulting coefficient 
diffusion, (extracted numerically from a linear fit of 
Fig.^Ja), can be shossai-.|to depend on the parameter K 
in the following wayE^oHj: 

lim (AJ 2 (t))/t oc K 5 /' 2 



V 

where t is the iteration time. 



(4) 




k 5/2 T 1/2 

FIG. 3. Diffusion rate for the discontinuous map (2) for 
different k and T versus k 5/2 T 1/2 . Open circles are for T = 1 
, full circles for k = 0.01 ; both for kT < 1. Full line is the 
best fit line D = 0.39fc 5/2 T 1/2 . 

In terms of the original variables one get D = 
(Al 2 (t))/t oc fc 5 / 2 VT (see Fig.|). The regime charac- 
terized by a diffusioiLjCoefficient scaling as K 5 / 2 has been 
called slow diffusionEl. This mark the difference from the 
standard quasilinear regime [K > 1) where-typically one 
has V ~ K 2 (superimposed to oscillationst-a) . 



The dependence on the power 5/2 has been explained^! 
in terms of a transport model based on a Markovian par- 
tition of the phase space for the STM and it clearly indi- 
cates that the random phase approximation!] cannot be 
applied in this case. 

Phases are indeed correlated within a time r (see Fig.|| 
b) during which the resulting motion cannot be chaotic 
nor diffusive. A correct evaluation of this time scale could 
be a key in understanding the "strange" exponent 5/2, 
even if numerically it is quite difficult to obtain sharp 
results. 



IV. QUANTUM MAP DYNAMICS 

Now we turn to the quantum analysis. Adopting the 
nowadays standard procedurea we analyze the quantum 
evolution of map (|^) by means of the one-period evolu- 
tion operator given explicitly by (h = 1): 



Ut = e 



-iTh 2 /2 -ik\ cos0| 



(5) 



where h = ~id/d6. Apart from the modulus in the 
potential V{6) — |cos0| the evolution operator Ut is 
exactly the same as the Kicked Rotator one. Anyway 
the presence of the modulus leads to many important 
physical differencies. Indeed, written in the momentum 
basis n, the matrix elements U n-m = (n\UT\m) decay as 



a power law \lA n 



l/\n 



away from the princi- 



pal diagonal (and not faster than exponentially as for the 
KRM). This case has been investigated for banded ran- 
dom matricesEj, where it was shown that cigenfunctions 
are also algebraically localized with the same exponent. 

The presence of power law localized eigenstates has 
major consequences. First of all it is not, a priori, obvious 
if the mechanism connected with the exponential dynam- 
ical localization holds even in this case. Moreover, while 
in case of exponential localization a unique measure of 
localization is defined (up to a constant), for algebraical 
localization different definitions of localizations can, in 
principle, give rise to different parametric dependences. 

Here we consider, as a measure of the degree of lo- 
calization, the variance £ CT of the stationary distribution 
P(n) = \il> n {t^ oc) " 



2. 



t = E^w] 1/2 (6) 

n 

which we expect to have a well defined classical limit. 

A. The Kicked Rotator Model for kT < 1 

Before we turn to numerical results, it is useful to con- 
sider, as a common reference, the KRM. This is far from 
being pedagogical, since the case we are interested in here 
(kT < 1) has not been an object of intense investigations 
in the past. Moreover, in this region, characterized by 
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classical regular motion, it is quite difficult to get homo- 
geneous results. For instance, one should expect that the 
localization length of the stationary distribution strongly 
depend on the initial conditions. Starting within a reg- 
ular region (stable island around periodic orbits) would 
result in a spreading width surely not larger than the size 
of the island (excluding the exponentially small tunnel- 
ing among different classical tori) . This was indeed what 
ShepelyanskyES found: the spreading of the quantum sta- 
tionary distribution can be roughly identified with the 
width ~ yfk/T of the main classical resonance, whose 
size is yk/T, when kT < 1. Another result was ob- 
tained again by Shepelyansky bv investigating directly 
the quasi-energy eigenfunctiongij : he found direct pro- 
portionality between their localization lengths and the 
parameter fc : I ~ fc/4. The different regimes in this 
undercritical case were also reported by IzrailevQ. Ac- 
cording to his Fig. 3 (see also text) the case K = kT < 1 
is marked by two different quantum borders. One is the 
condition for the applicability of common perturbative 
theory (k ~ 1) while the other one is the condition for 
the semiclassical approach to describe quasiperiodic or 
chaotic motion (Shuryak border k ~ T). When k < T 
the size of the nonlinear resonance is less than the dis- 
tance between neighboring unperturbed levels and the 
quasiperiodic classical behavior is suppressed by quan- 
tum effects. 

Our data on the behavior of £ CT confirm and extend this 
general picture. They can be summarized as follows: 

f) £ CT depends, for K = kT < 1, only on the scaling 
parameter k/T. 

2) as a function of k/T two different regimes can be 
numerically detected, one linear, when k < T (below the 
Shuryak border) and another, for k > T, where £ ff cx 

These results are presented in Figj4|. 

As one can see the scaling law is accurate up to eight 
orders of magnitude. The intersection point between the 
two lines ( k ~ T called Shuryak border), is the only 
transition point (for £ CT ) we are able to detect numeri- 
cally. Results have been checked to be independent from 
the choice of the initial state ip n = S n , no . Since typically 
the variance £ ff (i) = y (An 2 (t)) is an oscillatory func- 
tion of the iteration time t (quasi-periodic motion) the 
average (in time) value has been taken. Fluctuations are 
typically quite large, often on the same order of the av- 
erage value. On the basis of previous results we identify 
two different regimes, for the KRM in the classically reg- 
ular case kT < 1 : the perturbative one £ CT ~ k/T and 
the quasi-integrable one £ CT ~ y/k/T. 
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FIG. 4. Localization length as a function of the scaling pa- 
rameter k/T for the Kicked Rotator and kT < 1. Full symbols 
are obtained by fixing k and varying T : circles (k — 0.01), 
squares (k = 1), asterisks (k = 10). Open symbols are ob- 
tained by fixing T and varying k: circles (T = 0.01), squares 
(T — 0.1). Dashed line is 1.5fc/T (linear regime), while dot- 
ted line is 0.9 Jk/T. 



B. The discontinuous quantum model in the slow 
diffusive case 

In this section we analyze the quantum behavior of 
the discontinuous map The classical map is a good 
approximation to the real billiard dynamics only for 
e = kT/A -c 1. This means that we should consider, 
as a physical regime, only the slow diffusive one. The 
condition of applicability \l\/y/2E < 1 will be considered 
in details in section V. 

Iteration of the quantum map (||) in this regime typi- 
cally gives rise to the second moment (An 2 (t)) which is 
an oscillatory function of the time t. In particular, it is 
possible to characterize values of the parameters k and T 
which leads to periodic or irregular oscillations. In Fig.^ 
we show, for different k and T, the behavior of £ CT in time. 
As dashed lines we indicate the average values over few 
oscillation periods. We will back to this picture later on 
when we discuss the different quantum regimes. 

By varying k and T and taking £ CT as in FigJIl we ob- 
tain three different scaling regions: the first one for k < T 
in which £ CT ~ k/T, a second one for 1/VT > k > T 
characterized by ^ a ~ y/k/T and the third one, for 
1/VT < k < 1/T, in which & ~ D (see Fig.|). 
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FIG. 5. Wave packet spreading as a function of the it- 
eration time t. Here £ CT (t) = n2 |V'™(*)| 2 ] 1 ( a ) i s 
for k = 0.01 and T = 0.1 (perturbative region), dashed 
line is the average £ CT = 0.033; (b) is for k — 1 and 
T = 0.01 (quasi-integrable region), dashed line is the av- 
erage £ CT = 6.33; (c) is for k = 100 and T — 0.001 (dynamical 
localization region), dashed line is the average £ CT = 650. 

It is easy to identify the first two regions in close anal- 
ogy with those found for the KRM (see previous sub- 
section). While the existence of a perturbative region 
with the same characteristics of the KRM is not surpris- 
ing, much care must be taken in interpreting the quasi- 
integrable region. Indeed, while for the KRM we may 
properly speak about width of classical resonance, in this 
case we have no classical resonances at all. Nevertheless 
a close inspection of Fig.|| indicates the presence of is- 
lands of "quasi-integr ability" which means that trajec- 
tories spend a lot of time before leaving from them. This 
region is indeed dominated by classical Cantori, which 
act, from the quantum point of view as total barriers to 
the motion. These effects, namely the quantum propaga- 
tionJJorxjugh classical Cantori have been investigated in 
Rcf.El c3. In particular, a relation between width of the 
holes of Cantori and h should exist, in order to obtain a 



meaningful semiclassical limit. For instance Mackay and 
Meiss proposecS that Cantori could act as proper tori if 
the flux exchanged among different turnstiles is less than 

h. 
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FIG. 6. Localization length as a function of k/T. Three 
different sets of cases are shown, each one keeping T fixed and 
varying k: full circles (T = 0.001), open circles (T = 0.01) 
and asterisks (T = 0.05). Dashed line is fc/3T while dotted 
is y/k/T/2. Full lines are £ CT = D(T) for the three different 
sets, from the left to the right T = 0.05, 0.01, 0.001. 

Let us now analyze in more details the existence of a 
third scaling region for the localization length. The clas- 
sical analog of these quantum regions of localizations is 
a diffusive regime (after a transient time). Nevertheless 
the localization found is not "dynamical" in the sense 
that it cannot be derived for instance by taking the ap- 
proach used by Chirikov, Izrailcv and Shepelyansky for 
the KRMc-X Moreover, since this kind of localization is 
shared by the KRM, whose classical counterpart is reg- 
ular motion, it cannot be followed by classical diffusive 
excitation. A look at Fig.|| a-b confirms this view. 

On the other side, if the typical relation of the dynam- 
ical localization l a ~ D holds true, the diffusion rate, 
in order to produce an initial quantum classical-like dif- 
fusive spreading, has to be larger than the size of the 
classical "quasi-resonance" . This in turn allows us to 
estimate the second transition point as follows: 



D = D k' ,/2 VT = cyfkjT 



(7) 



which gives k cr — U c/DqT (where c and Dq are nu- 
merical constants of order one. One can then guess that, 
if any, the dynamical localization regime can exist for 
k > k cr = ^c/DqT. 
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FIG. 7. Stationary distribution averaged over few oscil- 
lation periods as a function of the momentum n. (a) is for 
k — 0.01 and T = 0.1 (perturbative region) (b) is for k = 1 
and T = 0.01 (quasi-integrable region) (c) is for k = 100 and 
T = 0.001 (dynamical localization region). Lines represent 
the power 1/n 4 and are shown to guide the eye. 

Our numerical computations indicate that the dynam- 
ical localization regime indeed exists, as can be inferred 
from Fig.[| The sharp rise of the full lines shown in FigJ|, 
(one for each T since the diffusion rate depends on T) in- 
dicates, without any doubt, the validity of the previous 
picture and the existence of the critical points, k ~ T 
and k ~ 1/VT. 

The existence of these thresholds for the Bunimovich 
stadium, and the regime of quasi-integrability as well, 
has been-confirmed analytically by Prange, Narevich and 
ZaitsevEl Indeed they were able to find an analytical 
expression for the eigenfunctions up to eyE ~ 1, which 
in terms of map variables reads k ~ 1/ \T. Above the 
threshold e\J~E ~ 1 the semiclassical perturbative ap- 
proach fails. According to our point of view this represent 
the k value necessary to start the classical-like diffusion 
process in presence of slow diffusion. 



Let us add few comments about the shape of the sta- 
tionary distribution. In all cases we have found good 
agreement with a power law distribution P(n) ~ \n — 
rto|~ 4 - In Fig.0 we show the correspondent stationary 
distributions for the cases of Fig.||. The lines, indicat- 
ing the power law behavior, are drawn to guide the eye. 
This means that results obtained from banded random 
matricestj can be extrapolated even when both random- 
ness and dynamical chaos are absent (perturbative and 
quasi-integrable regions). 

The presence of power law localized states for the dis- 
continuous map, also recently confirmed inEJ for the sta- 
dium eigenfunctions, should be somehow .put in relation 
with the exponential localization found inB for the eigen- 
functions of a rough billiard. This peculiarity should be 
in turn related with the particular boundary shape per- 
turbation. Indeed, while in the rough billiard considered 
irJa, a finite number M of harminics is necessary in order 
to produce the deformation from the circle, the stadium 
boundary perturbation is not analytical. The classical 
rough map is shownO to be chaotic when e > e c ~ M -5 ' 2 . 
It is then clear that when M — > oo KAM regular struc- 
tures disappear and the situation depicted by the dis- 
countinuous map appears. 

To end this section let us remark that a similar behav- 
ior happens if another definition of localization length is 
taken. In Fig.^J we show the inverse participation ratio of 
the quantum distribution ipr — 1/ ^ n |i/'n| 4 as a function 
of the scaling parameter k/T . The absence of a perturba- 
tive region is due to the fact that, by definition, ipr > 1. 
Moreover, let us note that a numerical constant has to 
be added in order to follow numerical data. Namely, the 
dynamical localization regime is marked by ipr ~ D/A. 
The simple proportionality between l a and ipr, even if 
not surprising (the same happens for the KRM), has to 
be considered accidental for power law localized distribu- 
tions. 
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FIG. 8. Inverse participation ratio ipr as a function of k/T. 
Three different sets of cases are shown, each one keeping T 
fixed and varying k: full circles (T = 0.001), open circles 
(T = 0.01) and asterisks (T = 0.05). Dotted line is O.Qy/k/T. 
Full lines are ipr = D(T)/4 for the three different sets, from 
the left to the right T = 0.05, 0.01, 0.001. 



V. BORDERS FOR THE STADIUM 

The results found in the previous sections can be ex- 
tended to the Bunimovich Stadium. This will lead to 
important estimates which enable us to discriminate be- 
tween different physical situations. In particular the crit- 
ical points found previously give rise to relations between 
the energy and the small parameter e. 




FIG. 9. Perturbative border for the Bunimovich Stadium. 
Nb is the quantum number corresponding to the eigen energy 
E. Solid circles represent numerical data, and straight line is 
the best-fit N b = 0.57/e 1 ' 04 . 



The first regime (perturbative) , is characterized by k < 
T or 2Ee < 1. From the billiard point of view this means 
that levels having an energy less than E ~ 1/e can be 
obtained perturbatively from those of the unperturbed 
spectrum (e = 0, e.g. the circle). The numerical study of 
the energy spectrum gives a confirmation for this border. 
We label the energy levels as Eiy(e). If e is sufficiently 
small we can assume the shift AE — -Ejv(e) — En(0) to 
be quite small. On the other hand the energy spectrum 
is characterized by an average level spacing SEn which 
is given approximately by the Thomas-Fermi formulate: 



SE 



N 



En 



iR 2 



(8) 



where to = 1 is the particle mass, R — 1 is the circle 
radius and Nj, is the number of levels up to energy E. In 
order to be in a perturbative regime there should be no 
levels overlapping, namely AEn(e) < 5En- This gives 
a relation between the energy E (or the level number 



Nb) and the small parameter e, which can be detected 
numerically. We show our numerical results in FigJ|. 
The best fit gives rise to 0.57/e 104 , which is in a good 
agreement with the theoretical prediction. Numerically, 
Ni, is obtained by comparing the eigenernergy list of the 
circular billiard and that of the perturbed Bunimovich 
stadium of e. _ . 

Moreover, a kind of semiclassical approach is possibles 
up to k = 1/Vr, or E ~ 1/e 4 (we omit a numerical con- 
stant in front of this expression : its value can be ob- 
tained only numerically). This means that the analytical 
approach to eigenvalues and eigenvectors is possible up 
to this energy value. Nevertheless, as soon as the map 
dynamics correctly approximates the real dynamics, the 
localization length should have a different energy depen- 
dence. Even in this case a direct numerical approach is 
needed in order to give a definite answer. 

Let us now come to the much more interesting case 
marked by the dynamical localization. We have found 
that this is possible only for E > 1/e 4 . According to the 
dynamical localization theory the quantum spreading will 
occur up to a time t B — D called break time. Following^ 
we may speak of a proper localized quantum regime only 
if is < terg, where t erg is the classical time in order 
to reach a stationary ergodic distribution, estimated as 
A; 2 ~ Dt B or t B ^ e~ 5 / 2 . Would t B be larger than 
t erg , classical and quantum distributions will both reach 
an ergodic stationary distribution. Putting t B = t erg or 
E ~ e -5 we get the last critical point above which we 
expect quantum as well classical ergodicity. 

The approach to ergodicity cannot be studied using 
this map. Indeed it represents a good approximation to 
real dynamics only for t <C t erg . This means that our 
results cannot be .compared with those found by Frahm 
and Shepelyanskynj about the approach to ergodicity via 
a Brcit-Wigner regime. This is a situation characterized 
by eigenfunctions delocalized on the energy shell but with 
many strong isolated peaks of probability. The presence 
of isolated peaks of probability in our case (see Fig. 7) is 
instead due to the classical phase space structure. 

While the ergodicity regime has been previously 
investigated^ using the NNLSD, the existence of two dif- 
ferent localized regimes (e _1 < E < e~ 4 and e~ 4 < E < 
e~ 5 ) has been only guessed on the basis of the similar- 
ity with the approximate map. In this case, the study 
of NNLSD should be probably accomplished by a direct 
study of eigenfunctions. This we argue, since it is not at 
all obvious how different kind of localization can affect 
the level statistics. A preliminary study in this direction 
can be found inEil. 

Let us note that, in order to correctly approximate real 
dynamics with the map, one has to require \J2E < 1, 
or IT < 2. One should then require that in both, quasi- 
integrable and dynamical localization regime, £ CT T < 2. 
This in turn means either 



= v IT 2 



(9) 
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(which is always justified since kT = 4e and e < 1) or 

fc 5/2 T 3/2 < 2 

which reduces to 

fc(fcT) 3 / 2 ~ £i? 1 /2 e 3/2 = E l/2 e 5/2 < 2 

or < 1 /e 5 which is the condition in order to have clas- 
sical ergodicity. This is the reason why we neglect this 
condition in Sec. IV. 

Most of the results were obtained by approximating 
the real quantum dynamics by means of a quantum map 
which is the quantum analog of a classical map approx- 
imating the classical dynamics. This kind of procedure 
is not new (see for instanceEH) . One may wonder if this 
is, at the end, close to the original model. The answer 
can come, of course, only from a direct numerical or ex- 
perimental, analysis of the quantum dynamics of wave 
packets inside the billiard. Nevertheless it is significa- 
tive that different approaches on the same modelEJ give 
results in good agreement with ours. 

After the completion of this wor.k|We became aware of 
other related works on the subjecttJ. In particular their 
numerical data, while confirming the existence of these 
regimes, indicate other different borders. However, from 
the map point of view, the transition at k = T (E ~ 1/e) 
is very sharp, we cannot exclude numerically the presence 
of a further border at k ~ T 1 / 3 (E ~ e~ 3 )- Indeed a close 
inspection at Fig. 6, indicates a very smooth transition 
toward the line £ CT — D. Further numerical calculations 
are required in order to show if quantum map also shows 
this border. 
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